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D Abstract 

Characteristics of the hadronic final state of diffractive deep inelastic scattering 
events, ep — > eXp, were studied in the kinematic range 4 < Mx < 35 GeV, 
4 < Q 2 < 150 GeV 2 , 70 < W < 250 GeV and 0.0003 < x P < 0.03 with the 
ZEUS detector at HERA using an integrated luminosity of 13.8 pb _1 . The events 
were tagged by identifying the diffractively scattered proton using the leading 
proton spectrometer. The properties of the hadronic final state, X, were studied 
in its center-of-mass frame using thrust, thrust angle, sphericity, energy flow, 
transverse energy flow and "seagull" distributions. As the invariant mass of the 
system increases, the final state becomes more collimated, more aligned and more 
asymmetric in the average transverse momentum with respect to the direction of 
the virtual photon. Comparisons of the properties of the hadronic final state with 
predictions from various Monte Carlo model generators suggest that the final state 
is dominated by qqg states at the parton level. 
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1 Introduction 

A class of deep inelastic scattering (DIS) events has been observed at HERA that have the 
characteristics of diffractive interactions. These events have a large rapidity gap M be- 
tween the recoil-proton system and the produced hadronic system, and a small momentum 
transfer to the proton ||. The events can be pictured in terms of the t-channel exchange 
of an object that carries the quantum numbers of the vacuum, called the Pomeron (IP), 
see Fig. |l](a). However, the nature of the Pomeron in DIS is at present far from clear. 
Measurements by the HI || £|, [5| and ZEUS f| collaborations have shown that, in QCD- 
inspired models of the diffractive process, the Pomeron can be described as an object 
whose partonic composition is dominated by gluons. Alternatively, the diffractive process 
can be described by the dissociation of the virtual photon into a qq or qqg final state that 
interacts with the proton by the exchange of a gluon ladder @. 

The study of the hadronic final state in e + e~ annihilation [Q has been a powerful tool 
in gaining information about the underlying partonic state. Similarly, the study of the 
partonic content of the hadronic final state in diffraction is a natural way to explore the 
dynamics of diffraction. 

In this paper, a study is reported of the hadronic system, X, produced in the DIS process 
ep — > eXp, where the diffractively scattered proton stays intact. The proton was de- 
tected and its three-momentum measured in the leading proton spectrometer (LPS) @. 
Diffractive events are defined, for the purpose of this paper, as those events which contain 
a proton with more than 97% of the initial proton beam energy. Previous results on 
hadronic final states in diffractive events at HERA have been obtained with the require- 
ment of a large rapidity gap between the observed hadronic system and the scattered 
proton |H| |llj. The results obtained with rapidity-gap events were either defined in a 
reduced phase space by imposing a cut in rapidity 0, [12| , or Monte Carlo simulated 



events were used to extrapolate the characteristics of diffractive events over the areas of 
phase space removed by the rapidity-gap cut [[Ll]]. By using the scattered proton to tag 
diffractive events, there is no need to rely on Monte Carlo generators to model correctly 
the part of the final state removed by the rapidity cuts, and the full angular coverage of 
the central detector can be used. 

The properties of the hadronic system, X, were studied in terms of global event-shape 
variables such as thrust and sphericity in the center-of-mass (CMS) frame of X. This 
is analogous to the studies of global event-shape variables in e + e~ annihilation || as a 
function of the CMS energy and to the analysis that led to the interpretation of three-jet 



events in terms of gluon bremsstrahlung [13]. In addition to global event-shape variables, 



the properties of the diffractive events were described in terms of inclusive distributions 
such as energy flow. 



2 Experimental setup 

The data were recorded in 1997 with the ZEUS detector and correspond to an integrated 
luminosity of 13.8 ± 0.3 pb _1 . A detailed description of the ZEUS detector can be found 
elsewhere [fL4 . A brief outline of the main detector components most relevant for this 



analysis is given below. Charged particles are tracked by the central tracking detector 
(CTD) ||15|| , which operates in a magnetic field of 1.43 T provided by a thin supercon- 
ducting coil. The CTD consists of 72 cylindrical drift chamber layers, organized in 9 
superlayers covering the polar angleQ region 15° < 9 < 164°. The transverse-momentum 
resolution for full-length tracks is <t(pt)/pt — 0.0058pr © 0.0065 © 0.0014/px, with p T 
in GeV. The high-resolution uranium-scintillator calorimeter (CAL) jl(J consists of three 
parts: the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each 
part is subdivided transversely into towers and longitudinally into one electromagnetic 
section (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test beam conditions, are a(E)/E = 0.18/ y/E for electrons 
and a(E)/E = 0.35/VE for hadrons (E in GeV). 

The LPS |J detects charged particles scattered at small angles and carrying a substantial 
fraction of the incoming proton momentum; these particles remain in the beampipe and 
their trajectory is measured by a system of position-sensitive silicon micro-strip detectors 
very close to the proton beam. The track deflection induced by the magnets in the proton 
beamline is used for the momentum analysis of the scattered proton. The LPS consists of 
six detector stations, SI to S6, placed along the beamline in the direction of the outgoing 
protons, at Z — 23.8 m, 40.3 m, 44.5 m, 63.0 m, 81.2 m and 90.0 m from the interaction 
point, respectively. In this analysis, only the stations S4 to S6 were used. These stations 
consist of two halves, each equipped with an assembly of six parallel planes of silicon 
micro-strip detectors which can be inserted to a position near the proton beam. Each 
detector plane has an elliptical cutout which follows the profile of the 10a envelope of 
the beam, where a is the standard deviation of the spatial distribution of the beam in 
the transverse plane. The acceptance of stations S4 to S6 of the LPS for protons close 
to the beam energy is a few per cent and approximately uniform for 0.075 < \t\ < 0.35 
GeV 2 . The LPS longitudinal momentum resolution is &(pz)/pz = 5 x 10~ 6 pz (pz in 
GeV). The transverse momentum resolutions, dominated by the proton beam emittance, 
are a Px = 35 MeV and a PY = 90 MeV. 



1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the center of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as r\ = — In (tan |), where the polar angle, 9, is measured with respect to the proton beam 
direction. 



3 Kinematic variables and event reconstruction 



3.1 Kinematic variables 

The event kinematics of DIS processes can be described by the negative squared four- 
momentum transfer at the lepton vertex, Q 2 = —q 2 = —(k — k) 2 (k and k denote 
the four- momenta of the initial- and final-state positrons, respectively), and the Bjorken 
scaling variable, x = Q 2 /(2P ■ q), where P is the four-momentum of the proton. The 
fraction of the energy transferred to the proton in its rest frame, y, is related to these two 
variables by y = P ■ q/(P ■ k) ~ Q 2 /xs, where y/s is the positron-proton CMS energy. 
The CMS energy of the j*P system, W, is given by W 2 = (q + P) 2 = Q 2 (l - x)/x + M 2 , 
where M p denotes the proton mass. Unless otherwise noted, the double angle method [17 
is used to reconstruct the kinematic variables from the measured quantities. 

Two more kinematic variables are needed to describe a diffractive DIS event, e(k) + 
p(P) — > e'(k') + X + p'(P'), where the scattered proton, p', has four-momentum P' = 
(E p i, P x , P Y , P' z ). The squared four- momentum transfer at the proton vertex is given by 

t = (P- P'f ~ -Pt, 

where P% = P'x + Py , and the fraction of the beam momentum retained by the final 
proton is 

\P'\ E p , 

L ~ \P\ ~ E ' 
I I v 

where E p = 820 GeV denotes the incident proton beam energy. Both t and xl are 

measured with the LPS. Other useful variables in diffractive DIS are: 

_ (P-P')-q M x + Q 2 -t 

XF ~ P-q ~ W 2 + Q 2 -M 2 ^ XL 

and 

Q 2 x Q 2 



P 



2{P-P')-q s P M x + Q 2 ' 

where M x is the invariant mass of the final-state, X, determined as described below. In 
the Pomeron-exchange picture, x-p is the fraction of the proton's momentum carried by 
the Pomeron. For a Pomeron with partonic constituents, /3 is then the fraction of the 
Pomeron's momentum carried by the parton that absorbs the virtual photon. 

3.2 Reconstruction of Mx 

Two methods were used in this analysis to determine the invariant mass of the system X. 
The first method measured the mass directly from X. The individual objects belonging 



to X were reconstructed by combining charged tracks measured in the CTD and neutral 
energy clusters measured in the calorimeter into energy flow objects (EFOs ||). In this 
procedure, the tracks and clusters associated with the scattered positron were not used. 
The four-vectors of the EFOs were calculated assuming the pion mass. The invariant 
mass can then be reconstructed using these four-momenta as 

m 2 x , efo = (£ E t f - (£ p Xi f - (£ p Yl f - (£ p Zl f. (i) 

i i i i 

The second method inferred Mx from the measurement of the final-state proton and the 
scattered positron, which enter into the reconstruction of the kinematic variables x and 
y via the double angle method. M x is then given by 

M 2 X , LPS = sy(l -x L -x). (2) 



Since the LPS method has better resolution at higher values of Mx and the EFO measure- 
ment is better at lower values, Mx was evaluated as the weighted average of the values 
obtained with the two methods. Before combining the two Mx values, correction factors 
obtained from Monte Carlo studies were applied. This resulted in a resolution on Mx of 
about 25% at low M x , improving to 15% at high M X - 



3.3 Boost into the 7*P center-of-mass system 

The boost into the CMS of X (7T), Fig. 0(b), was determined from the four-momenta 
of the 7* and P. The 7* four-momentum was calculated from the energy and angles of 
the scattered positron. The best reconstruction of the P four-momentum was obtained 
by combining the information from the LPS and the ZEUS main detector. The X and 
Y components of the IP four-momentum were taken to be equal to the negative of the 
X and Y momentum components of the scattered proton measured with the LPS. The 
Z component was taken to be the total Pz of all EFOs (P z ot = ^ Pzi) minus the Z 
component of the 7* four-momentum. Finally, the energy component was calculated by 
requiring that the 7*P invariant mass be equal to the value of Mx determined by the 
combination of the information from the LPS and the EFOs, as described earlier. 

The resolution on the angle between the 7*P axis and the ZEUS Z-axis, as estimated 
by Monte Carlo studies, was approximately inversely proportional to Mx and equal to 
10 mrad at Mx = 20 GeV. An improvement in the accuracy of the measurement by a 
factor of three was achieved using the LPS compared to a measurement without its use, 
i.e. substituting px = Py = 0. 



4 Event selection 

Deep inelastic scattering events were selected in a manner similar to those of the ZEUS 
F 2 analysis fL8[ by requiring: 



the presence of a fully contained and isolated positron with energy greater than 10 GeV 
and matched to a charged track, if within the acceptance of the CTD; 



a value of y greater than 0.05, evaluated using the Jacquet-Blondel method [IS]; 



a value of 8 = E tot - P^ ot in the range 40 < 5 < 65 GeV, where E tot and P tot 
are, respectively, the total energy and momentum of the event in the main detector, 
including the scattered positron; 

the Z coordinate of the event vertex within 50 cm of the nominal interaction point. 



Diffractive events were then selected by requiring a well reconstructed LPS track || 
carrying more than 95% of the beam energy (xl > 0.95). The LPS track was required 
to pass no closer than 0.04 cm to the HERA beampipe, and the quantity 5 + = E tot + 
P| ot + 2E p xl was required to be less than 2E p + 20 GeV to be consistent with longitudinal 
momentum conservation after taking the LPS resolution into account. Since the results 
presented here were found to be independent of t, no explicit cut on t was applied. 

To provide a sample for which the acceptance is large and uniform, to remove events 
measured with low resolution and to reject events corresponding to exclusive production 
of vector mesons, the following kinematic cuts were applied: 



70 < W < 250 GeV; 
0.0003 < a; P < 0.03; 
4 < Q 2 < 150 GeV 2 ; 
4 < M x < 35 GeV. 



In addition, at least four EFOs in the system X were required. These cuts define the 
kinematic range for all results presented in this paper, unless otherwise noted. The final 
data sample contained 2355 events. 



5 Models of the diffractive final state 

The data were compared with three Monte Carlo generators based on different theoretical 
models. For all generators, hadronization was simulated using the Lund string model as 
implemented in JETSET 7.4 [2(J. The models considered here produce either a qq pair 



or a qqg final state at the parton level. However, the dynamics of the production of these 
states is different for each model, thus yielding predictions that differ in their relative 
contributions of qq and qqg states, as well as in the final-state topology. Common to all 
three models is that the qqg-type events dominate the final state at high masses. In such 
events, the gluon usually travels in the direction of the Pomeron. 



5.1 Resolved Pomeron model 

In the resolved Pomeron approach pT| , the exchanged particle is assumed to have a par- 
tonic structure consisting of quarks and gluons. A sample of resolved Pomeron events was 
produced with the RAPGAP 2.08/06 generator [p2]] . The hadronic final state is simu- 
lated in analogy to ordinary DIS. The Pomeron parton density functions (PDFs) used were 
those determined by the HI collaboration from their measurement of d 3 a/dx-pdQ 2 d/3 |3|. 

When the virtual photon in a diffractive event interacts with a quark in the Pomeron, 
the lowest-order final state is a qq. At 0(a s ), qqg final states can be produced by the 
radiation of an extra gluon via the QCD-Compton (QCDC) process. In addition, qqg 
final states are produced via the boson-gluon fusion (BGF) process between the virtual 
photon and a gluon in the Pomeron. The Pomeron remnant consists of a quark in the qq 
and QCDC events and of a gluon in BGF events. Additional parton radiation from the 
remnant partons is suppressed. The following RAPGAP options were used: 

• the leading-order version of the HI Fit 2 was used. In the fit, more than 80 % of the 
Pomeron momentum is carried by the gluons. 

• the final-state parton system in gg-type events was simulated with limited transverse 
momentum with respect to the 7*IP axis (intrinsic transverse momentum, kx), ran- 
domly distributed according to exp(— 5.5/cf,) {kx in GeV); 

• to avoid divergences in the matrix elements for the 0(a s ) processes for massless quarks, 
the transverse momentum squared of any outgoing parton was required to be larger 
than 3 GeV 2 ; 

• higher-order QCD radiation was simulated using initial- and final-state parton showers 

(MEPS) @; 

6 



• the generator was run for light flavors and charm production (produced via the BGF 
process) separately, and the final sample was obtained by mixing the two according 
to their relative cross sections as predicted by RAPGAP. 

The curves corresponding to this model are labeled "RG resolved P" in later figures. The 
model labeled "RG qq part only" in later figures was obtained by selecting events from 
the above sample that were not produced via the BGF process. 



5.2 Photon dissociative model 

In the rest frame of the proton, diffractive scattering can be viewed as the dissociation of 
the virtual photon into a qq pair (Fig. |Tj(c)) well before the interaction with the proton. 
For higher masses, the qqg final state (Fig. |l](d)) becomes important. The dissociated 
photon system couples to the proton by color-singlet exchange. 

The simplest realization of this color singlet is the exchange of two gluons with opposite 
color charge (two-gluon exchange model). The RIDI2.0 program |24| implements a two- 
gluon model following the approach of Ryskin p5| , where the diffractive dissociation is 



treated in the framework of the leading log approximation (LLA) of perturbative QCD. 
The cross section is proportional to the square of the gluon density of the proton, which 



was taken from the CTEQ4M [26] parameterization of the proton structure function with 
an appropriately chosen cut-off parameter for the transverse momentum of the final-state 
gluon. The contributions of both transversely and longitudinally polarized photons are 
included. A large theoretical uncertainty in the relative contributions of qq and qqg states 
remains. This uncertainty is reflected in the wide range allowed for the (a^-dependent) 
fc-factors, which determine the relative cross sections of the qq and qqg contributions. 



m 



Another model based on photon dissociation is that of Golec-Biernat and Wusthoff p?| 
which the virtual photon splits into a qq or qqg color dipole. The interaction of this dipole 
with the proton can be described by an effective dipole cross section taking the dynamics 
of saturation into account. The dipole cross section was parameterized using the HERA 
measurements of the DIS total cross section between Q 2 = 0.1 GeV 2 and Q 2 = 100 GeV 2 . 
The production of qq and qqg final states is calculated using these parameters. The 
relative contribution of qq and qqg events is fixed by QCD color factors. This model 
is implemented in the SATRAP Monte Carlo generator |2£| which uses the additional 



approximation of strong ordering in the transverse momenta: kx{g) "C kx{q, q)- In the 
course of the ZEUS analysis of diffractive 3-jet events |29|1 , it was found that the modeling 
of higher-order QCD processes was inadequately implemented in SATRAP. Most notably, 
no initial-state parton cascades were included, and the final-state QCD radiation from 
the gluon in the dominant qqg contribution was suppressed. A new implementation of 



higher-order QCD processes in SATRAP was carried out [30, 31 1, in which the color-dipole 
model (CDM) was implemented in a similar fashion to that in RAPGAP. This model is 
referred to as SATRAP-CDM. 



5.3 JETSET 

If the diffractive DIS hadronic system, X, is produced by a virtual photon splitting into 
a qq pair, its properties would be expected to be similar to those of the hadronic system 
produced by e + e~ annihilation at a CMS energy y^ = Mx- The e + e~ final state, as 
simulated by the JETSET program [pOR , was used as a baseline to which both the data 
and the diffractive Monte Carlo generators were compared. JETSET is known to describe 
accurately many details of the final state in e + e~ collisions. 



6 Data correction and background 

Monte Carlo simulations were used to correct the data for the resolution and acceptance 
of the main detector and the LPS, and to estimate the size and influence of the back- 
ground. To estimate the model dependence of these corrections, two different generators 
were used: RAPGAP [F2^] and RIDI @- RAPGAP events were generated with the HI 



QCD Fit 2 || for the IP structure function, as described in Section |5.1| , except that 
CDM rather then MEPS was used for higher-order initial- and final-state QCD radia- 
tion. This model gives a poor description of the hadronic final state as well as of some 
kinematic variables as measured at the detector level. To obtain a sample which gives a 
good description of all observed distributions, the RAPGAP generator was reweighted to 
reproduce the data. The corresponding distributions are labeled "mod. RG" in Figs. to 
|^. A similar reweighting procedure was carried out for RIDI. Both samples were passed 
through a GEANT 3.13 |32| simulation of the ZEUS detector, subjected to the same 
trigger requirements as the data and processed by the same reconstruction programs. 

In addition, changes in the HERA beamline parameters during the running period and 
the finite resolution of the motors that determine the position of the LPS, neither of which 
was simulated in the Monte Carlo, were taken into account by reweighting and smearing 
the simulated xl distribution. 

The primary source of background in this analysis is the accidental overlap of a DIS event 
with an unrelated beam-halo proton measured in the LPS. This contribution was assumed 
to be constant over the running period. To estimate its size, a sample of such background 
events was studied. Background protons were identified by selecting DIS events with 
E tot + P|°* > 100 GeV (using only CAL) and an LPS track with xl > 0.9. From energy 
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conservation allowing for detector resolution, such events must result from an accidental 
overlap. The LPS information from these events was then combined with a sample of 
non-diffractive DIS Monte Carlo events generated using DJANGOH |£|. The resulting 



sample was normalized to match the upper tail of the 5 + distribution for the data shown 
in Fig. 0. With this method, the level of background contamination, after all selection 
cuts, was estimated to be 5.2%. Its contribution was statistically subtracted in all results 
presented below. 

Figure |3| shows the distributions of Q 2 , W, xl and —t for events which pass all selection 
cuts except the ones imposed on the variables shown (indicated by the arrows). Also 
shown in the figure are the distributions from the modified RAPGAP simulation plus 
background which was normalized so that the weighted sum of the modified RAPGAP 
and the background Monte Carlo events describes the Q 2 distribution in the selected 
range, see Fig. ||| (a). The kinematic distributions of the data, with the exception of 
the t distribution, are well described by the modified RAPGAP Monte Carlo sample in 
combination with the background sample. 

Figure f| (a-c) shows some of the properties of the hadronic final state: the distribution 
of the number of EFOs, the EFO energy in the CMS frame of the 7* IP system and r] max , 
where r] max is the pseudorapidity of the most-forward calorimeter cluster of energy greater 
than 400 MeV. The data are well described by the modified RAPGAP sample in combina- 
tion with the background sample. Figure ^ (d) shows the average Mx reconstructed using 
the LPS (Eq. (^J)) as a function of the invariant mass calculated using EFOs (Eq. ([I])). 
The correlation seen in the data is well reproduced by the Monte Carlo simulation. The 
fact that the correlation is approximately linear in the selected mass range confirms that 
the final state is well contained in the ZEUS detector. 

Figure || shows the distributions of M x , ip and (3. After the reweighting procedure, 
all data distributions are well described by the modified RAPGAP Monte Carlo sample. 
This reweighting, which mainly affected the low Mx and high ip regions, had little effect 
on the /3 distribution which, although peaking at low values, is well described by the 
relatively flat Pomeron PDFs used in RAPGAP. RIDI (not shown in Figs. H to |5|) also 
reproduces the data satisfactorily. 



7 Systematic uncertainties 

The systematic uncertainties were obtained by studying the sensitivity of the results to 
the uncertainties in the understanding of the ZEUS detector response and by modifying 
the analysis procedures as listed below. 

Uncertainties related to the understanding of the ZEUS main detector include: 



• the uncertainty on the absolute energy scales in each major section of the calorimeter 
(FCAL, BCAL, and RCAL), which for this analysis was understood to a level of 3% 
for hadrons and 2% for the scattered positron; 

• the uncertainty on the survey measurements of the positions of the major sections 
of the calorimeter with respect to each other and to the HERA beamline, which are 
accurate to 1-2 mm. 

Uncertainties related to the LPS include: 



• possible shifts in xl, Px, and py of ±0.003, ±10 MeV, and ±50 MeV, respectively as 
determined from an analysis of elastic p photoproduction (ep — > epp); 

• the uncertainty on the beampipe position, as determined by alignment studies, was 
taken into account by changing the cut on the distance of closest approach of the 
proton track to the beampipe by ±400 //m; 

• the uncertainty in the Monte Carlo simulation of the proton track reconstruction was 
taken into account by applying tighter cuts on the quality of the Monte Carlo track. 

To check for the effect of possible particle losses into the forward beampipe, and to 
account for the differences in modeling this forward region in the Monte Carlo generators, 
the energy deposited in the inner ring of the FCAL, which covers approximately the 
pseudorapidity range of 2.7 < r\ < 3.9, was scaled by ±25%, as suggested by MC studies. 

The normalization of the background from an accidental overlap of a DIS event with an 
unrelated proton was estimated using the 5 + distribution of Fig. |^. The uncertainty caused 
by the background subtraction was determined by changing the background normalization 

by ±50%. 

To estimate the model dependence of the results, the data were also corrected using RIDI. 
The difference between the RIDI and modified RAPGAP results was taken as an estimate 
of the model uncertainty. This uncertainty was assumed to be symmetric with respect to 
the nominal results obtained using the modified RAPGAP sample. 

The largest contributions to the systematic uncertainty typically originate from the model 
dependence and the uncertainty in the xl reconstruction. All systematic uncertainties 
were assumed to be independent and were calculated separately for positive and negative 
variations with respect to the nominal value. The total positive and negative systematic 
uncertainties were calculated as the corresponding sums in quadrature. 
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8 Global event-shape variables 

The results presented here are corrected to the hadron level in the kinematic range defined 
in Section |] with at least four stable particles in the final state. A particle is considered 
stable if its lifetime is longer than 3 x 10~ 10 s. If the lifetime is shorter than 3 x 10~ 10 s, 
the daughter particles (with lifetime > 3 x 10~ 10 s) are included. This corresponds to the 
final-state definition used by e + e~ experiments 0. 

8.1 Thrust and sphericity 

The event-shape- variable thrust (T) is calculated by determining the unit vector n which 
maximizes 

T(A): £<|a-wl 



£» \Pi 



where Pi is the three-momentum of a final-state particle and the sum is over all particles 
belonging to the system under study. The resulting axis h is called the thrust axis and 
T{p) is the thrust value. For collimated two-jet events, the value of T approaches 1, while 
events with an isotropic shape yield values close to 0.5. 



The sphericity (S) is defined as 



S = ^(A 2 + A 3 ), 



where A 2 and A3 are the two smallest eigenvalues (corresponding to orthogonal eigenvec- 
tors) of the sphericity tensor 

S^ = £if^ a,(3 = x,y,z. (3) 

For isotropic events, S approaches 1, and for collimated two-jet events, S is close to 0. 

The sphericity and thrust distributions in the CMS of the final state in diffractive DIS 
are shown in Fig. ^| for two bins of Mx, and are compared to the ones observed in e + e~ 
annihilation from the TASSO collaboration J8| at y/s —(Mx)- The diffractive events show 
thrust (sphericity) distributions which are broader and shifted to lower (higher) values 
compared to the e + e~ results, indicating that they are more isotropic. The diffractive 
events become more collimated as Mx increases, a trend also observed in e + e~ annihila- 
tion. If the virtual photon in diffractive DIS fluctuates only into a qq state, the resulting 
hadronic final state should develop in a manner similar to e + e~ — > qq at a CMS energy 
a/s —(Mx)- Deviations from this e + e~-type of behavior are expected, however, since 
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there is a significant contribution from the qqg diagram in the color field of the proton 
(Fig. 0(d)), which is not present in e + e~ collisions, where gluons can only be produced 
via higher-order QCD radiation. 

Figures (a) and (c) show the average values of thrust, (T), and sphericity, (S), measured 
in six bins of Mj. The value of (T) increases and (S) decreases with increasing M% with 
slopes similar to those found in the e + e~ data ||. On average, the diffractive events are 
less collimated than e + e~ events at a similar CMS energy. 

Figures |7| (b) and (d) show the same data, but now compared to several Monte Carlo 
modelsQ. Also plotted in Fig. [7] (b) are the average thrust values measured by the HI 



collaboration using a diffractive DIS sample tagged with rapidity gaps fTT |. The present 
measurements are in good agreement with the HI results which have been corrected to the 
full phase space by a MC simulation. The results are also consistent with those presented 
in a previous ZEUS publication [|T0H , which were obtained for a small range in rj. 



The Monte Carlo models of diffractive scattering reproduce the general trend of the data. 
RAPGAP and SATRAP-CDM give a good description of (T) over almost the full M x 
range, while RIDI fails for Mx < 20 GeV. The diffractive models produce events more 
isotropic than the ones generated with JETSET and measured in e + e~ annihilation. This 
is mainly due to the inclusion of qqg final states. Figures [7| (b) and (d) also show the 
RAPGAP contribution for the qq diagram alone. These curves demonstrate the need for 
the qqg contribution, especially at the higher Mx values. 

The final state measured in terms of thrust and sphericity has no strong dependence on any 
kinematic variable other than Mx- Figure [8] shows (T) as a function of Mx for xp < 0.01, 
where Pomeron exchange dominates, and for xp > 0.01, where Reggeon exchange may 
become more important. No significant differences are observed and all models give a 
reasonable description of the data. The data were also split into two samples of Q 2 , t and 
x; the dependence on these variables was less pronounced than that for x-p. 



8.2 Transverse momenta in and out of the event plane 

The shape of the system X was also studied in terms of two more variables, Pj~ in and 
Pt out) w hich measure the transverse momentum in and out of the event plane. The event 
plane is defined by the eigenvectors of the sphericity tensor associated with the two largest 
eigenvalues, Ai and A2 (see Eq. (|])). These transverse momenta have also been studied in 
e + er experiments, since they are sensitive to gluon bremsstrahlung. They are evaluated 



2 Henceforth, all comparisons of the data are made to MC models that have not been reweighted in the 
manner described in Sec. g. 
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by multiplying the average momentum squared of the n particles in the event by the two 
smallest sphericity eigenvalues A2 and A3, respectively: 

Phn = ^A 2 , (4) 

Pr,out — — fr^ 3 ' ( 5 ) 

where pj is the squared momentum of the zth particle in the event. By definition, P^ in > 
Pj-out- The observation of a difference in the Mx behavior of these two variables could 
be explained by the presence of events of planar shape, such as events with three partons 
in the final state. 

Figures |] (a) and (b) show the average values of P% in and P% out for the diffractive DIS 
sample as a function of M X - On average, P| in is about a factor of three larger than P% out . 
The dependence of P% in on Mx is almost linear over the range studied, while the Mx de- 
pendence of P| out becomes weaker as M x increases, indicating that the diffractive events 
become more planar. Figure |9] also shows the predictions of the diffractive DIS genera- 
tors and JETSET. The RAPGAP and SATRAP-CDM samples give a good description of 
P^in for low masses, but underestimate its value for M x ^ 15 GeV. RIDI overestimates 
P% in for most of the measured Mx region. Transverse momentum out of the event plane 
can be produced by the fragmentation process. While RAPGAP and SATRAP-CDM 
give a reasonable description of the data, RIDI slightly overestimates the data at low Mx 
values. The transverse momentum distributions in and out of the event plane predicted 
by JETSET are in excellent agreement with the diffractive scattering data. 

In summary, the diffractive hadronic final state shows trends in (T) and (S) different 
from those observed in e + e~ annihilation. These differences can be explained by the 
presence of qqg events from BGF (RAPGAP) or from photon dissociation (SATRAP- 
CDM, RIDI) which are absent in e + e~ annihilation and JETSET, where the only source 
of gluons in the final state is from QCD radiation. The transverse momentum out of the 
event plane obtained from diffractive Monte Carlo models and JETSET reproduces the 
data, indicating that the hadronization process is consistently modeled and is similar in 
diffractive DIS and e + e~ annihilation. 



9 Thrust axis orientation 



The orientation of the entire hadronic final state, X, with respect to the 7*IP axis can be 
studied by measuring the angle thrust (as shown in Fig. [J] (b)) between this axis and the 
thrust axis. 
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The values of ^Thrust and the transverse momentum, P^ hrust , relative to the thrust axis, 
defined by 

pThrust^^g^^ (6) 

are sensitive to various sources of transverse momentum. These include the intrin- 
sic transverse momentum, kr, of the partons in the proton, which was found to be 
1.69±0. 18^20 GeV in a recent ZEUS publication pl |, intrinsic transverse momentum 



of the partons in the Pomeron, and contributions of hard processes, e.g. final-state gluon 
emission [B31. 



Figure [U] (a) shows the average thrust angle (6 Thl . ust } as a function of M X - The angle is 
largest at small masses and falls off steeply as Mx increases, indicating that the event 
becomes more aligned with the 7* IP axis as the CMS energy increases. This general trend 
is described by all Monte Carlo models shown, but only RIDI produces a thrust angle as 
large as that observed in the data. SATRAP-CDM and RAPGAP fall about equally short 
of the observed angles. Given that the hadronization is well described by all Monte Carlo 



models (see Section |3.2|) , the measurement of (^Thrust) implies that additional sources 
of transverse momentum as discussed above must be important. Indeed, the approach 
implemented in RIDI favors the production of partons with relatively large transverse 



momenta, typically of the order of 1 GeV |24| . 

Figure |Il] (b) shows that the average transverse momentum, (p^ hrust ), produced in diffrac- 
tive scattering is almost independent of Mx for Mx > 10 GeV, with an average value of 
about 2 GeV. In this region of Mx, the independence is correctly reproduced by all Monte 
Carlo models, but the RAPGAP prediction is too low. The need for the qqg contribution is 
illustrated by the curve showing only that part of the RAPGAP prediction corresponding 
to the qq final state. 



10 Energy flow 

Another measurement of the event topology is the distribution of energy as a function 
of the pseudorapidity of a hadron with respect to the 7* IP CMS axis. This distribution, 



commonly referred to as the energy flow, is shown in Fig. |TT] for data in three Mx ranges. 
The data for low Mx have a Gaussian-like shape, but for Mx > 7.5 GeV a structure with 
two peaks develops. This structure becomes more pronounced as the mass increases. A 
slight asymmetry develops in the data, with more energy being produced in the 7* hemi- 
sphere. For comparison, the predictions from the RAPGAP, SATRAP-CDM and RIDI 
Monte Carlo generators are also shown. The SATRAP-CDM and RAPGAP generators 
predict too broad a rapidity distribution, and display a separation between the 7* and 
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P hemispheres in the lowest Mx bin that is not exhibited by the data. RIDI gives a rea- 
sonable description of the data in the lowest Mx bin but has a different shape at higher 
Mx- The asymmetry indicated in the data is slightly larger than that produced by the 
diffractive Monte Carlo models. For those values of \q\ at which the HI collaboration has 
also published data ||, the energy flow is in good agreement, except at the highest M x 
where the HI data are somewhat narrower. 

Figure |TJ shows the distribution of transverse energy, E T = ^2 Ei sin 0j, as a function of 
rj. Similar discrepancies between the data and the Monte Carlo events were observed for 
the transverse-energy flow as was observed for the energy-flow distribution. 



11 Seagull distribution 

The small asymmetry between positive and negative rj observed in the energy-flow plot 
of Fig. |ll| can be further investigated using the distribution of transverse momentum of 
the particles belonging to the system X. 

In inclusive DIS, ep — > eX, the fragmentation of the hadronic system occurs between the 
struck quark (forming the system X) and the remaining quarks in the proton (forming 
the proton remnant). Measurements in the Breit frame [36| have shown that, whereas the 



particle multiplicity and momentum distributions in the hemisphere of the struck quark 
are roughly consistent with those measured in e + e~ — ► qq, particles are produced with 
smaller average transverse momentum in the proton-remnant hemisphere. 

In diffractive DIS, an asymmetry in the momentum distribution between the 7* and 
IP hemispheres could be observed if the Pomeron behaves as an extended object like 
the proton and produces a remnant after the scattering process. This asymmetry is 
usually investigated by studying single-particle distributions as a function of the scaled 
longitudinal momentum, xf, and the momentum transverse to the photon direction, px- 
The photon direction defines the longitudinal axis in the 7*IP CMS as well as in the 7*p 
CMS. If pt and p^ are the momentum components of a final-state hadron perpendicular 
and parallel, respectively, to this axis, the variable xf is given by: 



max 



xf = Pl/p 

where positive xf is in the direction of the 7* and p™ ax is the maximum kinematically 
allowed value of p L . In the 7*IP CMS, p^ ax = M x /2. For the YP case, p^ ax = W/2. 

The average p T of particles as a function of xp, commonly referred to as the "seagull 
plot", is shown in Fig. [13] for 11 < Mx < 17.8 GeV. Also plotted are 7*p data from the 
EMC [37] collaboration at W — 14 GeV, equal to the average value of Mx in this bin. 

15 



The EMC DIS up — > fiX data indicate a suppression of the average p\ associated with a 
proton remnant which is not as apparent in the diffractive data in this mass range. 

Shown in Fig. [T3| (a)-(c) is the seagull plot for three different Mx bins compared with 
the predictions from RAPGAP, SATRAP-CDM and RIDI. The data exhibit a growing 
asymmetry as Mx increases. This asymmetry can also be seen in the ratios of the av- 
erage squared transverse-momentum in the 7* and IP hemispheres as a function of \xp\ 
(Fig. H (d)-(f)). The data are well reproduced by both RAPGAP and RIDI, while for 
Mx > 16 GeV, the transverse momentum generated by SATRAP-CDM is smaller than 
in the data. RAPGAP and SATRAP-CDM describe the size of the asymmetry correctly 
in all Mx bins, while for Mx > 16 GeV RIDI slightly underestimates the transverse 
momentum in the IP direction, resulting in a slightly larger asymmetry than that seen in 
the data. RAPGAP produces this asymmetry by including a Pomeron remnant. RIDI 
and SATRAP-CDM, in contrast, produce the asymmetry via the qqg diagram. 



12 Conclusion 

A study of the hadronic system, X, in the reaction ep — > eXp has been reported for 
the kinematic range 4 < M x < 35 GeV, 4 < Q 2 < 150 GeV 2 , 0.0003 < x F < 0.03 
and 70 < W < 250 GeV. The use of the LPS allows diffractive events to be tagged 
without applying cuts on the system X. It also provides a powerful constraint on the 
diffractive kinematics, allowing, for example, an accurate determination of the 7*P axis 
in the center-of-mass frame of the system X. 

The diffractive hadronic final state becomes more collimated as the invariant mass Mx 
of the system increases. This trend is similar to the one observed in e + e~ annihilation. 
However, on average the diffractive final state is more isotropic. This can be attributed 
to contributions not present in e + e~ annihilation, such as the boson-gluon fusion process 
in the resolved Pomeron approach, or qqg production from the dissociation of the virtual 
photon. 

The mean transverse momentum out of the event plane is similar to that found in e + e~ 
annihilation, indicating the universality of the hadronization. Even after considering the 
broadening effects of hadronization, it is apparent that more kx than is usually associated 
with the resolved Pomeron Monte Carlo approach is required to accommodate the large 
thrust angle and narrow energy flows at low Mx- 

Particle production becomes asymmetric along the 7*P axis as Mx increases, resulting in 
more average transverse momentum in the virtual-photon hemisphere. This asymmetry 
is consistent with both the concept of a remnant in the resolved Pomeron model and with 
the production of qqg final states in the photon-dissociation approach. The comparison 
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with the Monte Carlo models suggests that a dominant gluon contribution to the partonic 
final state is necessary. 

The invariant mass, Mx, of the hadronic system is the only variable upon which the 
characteristics of the system was found to depend. The system is independent of the 
momentum transfer at either the positron vertex, Q 2 , or at the proton vertex, t. Neither 
is there any dependence on x or on the fractional momentum of the Pomeron, x-p. 

Many models of diffraction are able to reproduce the measured diffractive cross sections. 
However, none of the models discussed here is able to describe all aspects of the data. It 
is clear, therefore, that measurements of the detailed characteristics of diffractive events, 
such as presented in this paper, will become more and more crucial in understanding the 
underlying physics of diffraction in deep inelastic scattering. 
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Figure 1: Different representations of diffractive DIS: (a) basic Pomeron-induced 
picture of diffraction; (b) the same process, in the 7*IP center- of-mass frame. The 
thrust angle, 9 = ^Thrust; is defined as the angle between the event axis and the 7*IP 
axis; (c) diffraction viewed in the rest frame of the proton as the fluctuation of the 
virtual photon well before the interaction with the proton to a qq pair; and (d) into 
a qqg system. 
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Figure 2: Normalization of the background. The quantity 5 + = T,(Ei + Pzi) + 
2E p Xl is shown for data (points, with statistical error only), background (shaded 
histogram) and modified RAP GAP plus background (solid histogram). The normal- 
ization was obtained as described in the text. The vertical arrow indicates the cut 
below which events were accepted for the analysis. 
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Figure 3: A comparison of data and Monte Carlo simulation for distributions of 
the kinematic variables of the diffractive DIS sample: (a) Q 2 , (b) W , (c) xl, and 
(d) —t. The data are displayed as points, with statistical errors only, and are com- 
pared to the modified RAPGAP plus background (solid histogram) and background 
alone (shaded histogram) . The normalization was obtained as described in the text. 
Vertical arrows in (a) and (b) indicate the values ofQ 2 and W between which events 
were selected for this analysis. Events with xl values below that indicated by the 
arrow in (c) were rejected. No cut on t was imposed. 
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Figure 4: A comparison of data and Monte Carlo simulation for distributions 
related to the properties of the reconstructed hadronic final state: (a) the number of 
EFOs; (b) the energy spectrum of the EFOs in the 7*IP frame; (c) i] max ; and (d) the 
average invariant-mass Mx measured with the LPS alone versus Mx determined 
from the EFOs. The data are displayed as points, with statistical errors only, 
and are compared to the modified RAPGAP plus background (solid histogram), 
background alone (shaded histogram) or modified RAPGAP alone (curve in (d)). 
The vertical arrow in (a) indicates the cut above which events were selected for this 
analysis. The region of Mx considered is shown by the arrows in (d). No cuts were 
imposed on the particle energy and r) max . 
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Figure 5: A comparison of data and Monte Carlo simulation for measured 
distributions of (a) Mx, (b) ip and (c) [3. The data are displayed as points, with 
statistical errors only, and are compared to the modified RAPGAP plus background 
(solid histogram). The background estimate is shown by the shaded histogram. The 
vertical arrows in (a) and (b) indicate the values of Mx and x-p between which 
events were selected for this analysis. No cut on (3 was imposed. 
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Figure 6: Thrust, T, and sphericity, S, distributions of the diffractive DIS 
hadronic final state (filled circles) compared to distributions measured in e + e~ col- 
lisions by the TASSO collaboration jffl (open circles) at the same CMS energies of 
the system X. In (a) and (c), results for 11 < M x < 17.8 GeV ((M x )= 14.0 GeV) 
are compared to measurements in e + e~ collisions at ^/s = 14 GeV; in (b) and (d), 
the results for 17.8 < M x < 27.7 GeV ((M x )= 22.0 GeV) are compared to e+e" 
results at s/s = 22 GeV. The inner error bars show statistical uncertainties only; 
the outer bars show the statistical and systematic uncertainties added in quadrature. 
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Figure 7: Average thrust (T) and sphericity (S) of the diffractive DIS hadronic 
final state as a function of Mx- Shown for comparison are results from e + e~ 
collisions and predictions from three Monte Carlo generators, RAPGAP (solid), 
RAPGAP qq only (dot-dashed), SATRAP-CDM (dashed) and RIDI (dotted). Also 
shown in (b) are thrust values from the HI Collaboration for diffractive DIS events 
tagged using rapidity gaps / (Il|/ . The inner error bars show statistical uncertain- 
ties only; the outer bars are the statistical and systematic uncertainties added in 
quadrature. 
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Figure 8: (a) Average thrust, (T) , of the diffractive DIS hadronic final state as 
a function of Mx for two different samples of events: x-p < 0.01 (solid circles) and 
xp > 0.01 (open circles), (b) Ratio of average thrust values T< = T(xp < 0.01) 
and T> = T(x$> > 0.01) as a function of M x for data (solid circles), and for 
predictions from RAPGAP (solid curve), SATRAP-CDM (dashed curve), and 
RIDI (dotted curve). The inner error bars show statistical uncertainties only; the 
outer bars are the statistical and systematic uncertainties added in quadrature. 
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Figure 9: (a) Average P^ in and (b) P^ out of the diffractive DIS hadronic final 
state as a function of Mx- Shown for comparison are predictions from three diffrac- 
tive Monte Carlo generators, RAPGAP (solid curve), SATRAP-CDM (dashed 
curve), and RIDI (dotted curve) together with e + e~ results as simulated by JETSET 
(dashed- dotted curve). The inner error bars show statistical uncertainties only; the 
outer bars are the statistical and systematic uncertainties added in quadrature. 
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Figure 10: (a) The average value of the angle ^Thrust of the thrust axis in the 
7*IP frame and (b) the corresponding transverse momentum p^ hrust compared to the 
prediction of RAPGAP (solid curve), qq part of RAPGAP (dot-dashed), SATRAP- 
CDM (dashed curve) and RIDI (dotted curve). The inner error bars show statistical 
uncertainties only; the outer bars show the statistical and systematic uncertainties 
added in quadrature. 
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Figure 11: The energy flow as a function of pseudorapidity, i], in the 7*IP- 
CMS, for various ranges of Mx ■ The data are shown as solid points, while the 
Monte Carlo predictions from RAPGAP, SATRAP-CDM and RIDI are shown as 
solid, dashed and dotted curves, respectively. Positive r\ is in the direction of the 
7*. The inner error bars show statistical uncertainties only; the outer bars show 
the statistical and systematic uncertainties added in quadrature. Also shown are 
the results from the HI collaboration (open squares) obtained from a measurement 
based on large rapidity gaps. 
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Figure 12: The transverse energy flow as a function of pseudorapidity, rj, in the 
^IP-CMS, for various ranges of Mj. 27ie data are shown as solid points, while 
the Monte Carlo predictions from RAPGAP, SATRAP-CDM and RIDI are shown 
as solid, dashed and dotted curves, respectively. Positive rj is in the direction of the 
7*. The inner error bars show statistical uncertainties only; the outer bars show 
the statistical and systematic uncertainties added in quadrature. 
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Figure 13: Average squared transverse momentum of particles measured in the 
center- of -mass frame of the system X as a function of xp for diffractive events 
from this analysis (solid circles) with 11 < M x < 17.8 GeV ({M x )= 14.0 GeV). 
Also shown (open circles) is the same quantity for inclusive DIS fip — > \xX data 
from the EMC collaboration fislj] at W = (Mx)- Positive xp is in the direction of 
the virtual photon. 
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Figure 14: Average squared transverse momentum as a function of xf (seagull 
plot) in three bins of Mx in the ^IP-CMS (a-c) and the ratios of average momenta 
in the 7* and IP hemisphere as function of \xp\ (d-f). The data (points) are com- 
pared to three models: RAPGAP (solid curve), SATRAP-CDM (dashed curve), 
and RIDI (dotted curve). Positive xf is in the direction of the virtual photon. 
The inner error bars show statistical uncertainties only; the outer bars show the 
statistical and systematic uncertainties added in quadrature. 
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